Matter-wave 2D solitons in crossed linear and nonlinear optical lattices 
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It is demonstrated the existence of multidimensional matter-wave solitons in a crossed optical 
lattice (OL) with linear OL in the x— direction and nonlinear OL (NOL) in the y— direction, where 
the NOL can be generated by a periodic spatial modulation of the scattering length using an 
optically induced Feshbach resonance. In particular, we show that such crossed linear and nonlinear 
OL allows to stabilize two-dimensional (2D) solitons against decay or collapse for both attractive and 
repulsive interactions. The solutions for the soliton stability are investigated analytically, by using 
a multi-Gaussian variational approach (VA), with the Vakhitov-Kolokolov (VK) necessary criterion 
for stability; and numerically, by using the relaxation method and direct numerical time integrations 
of the Gross-Pitaevskii equation (GPE). Very good agreement of the results corresponding to both 
treatments is observed. 

PACS numbers: 67.85.Hj, 03.75.Lm, 03.75.Kk, 67.85. Jk 



INTRODUCTION 

Bose-Einstein condensates (BEC) in optical lattices 
(OL) are presently attracting a great deal of interest, 
both experimentally and theoretically [3, 0] , in connec- 
tion with a series of physical phenomena which arise in 
condensed matter physics, these including Bloch oscilla- 
tions 3-5], generation of coherent atomic pulses (atom 
laser) |6|K dynamical localization 0-111, Landau- Zener 
tunneling (Iol - [l2l |. superfluid Mott transitions [l3|, etc. 
The flexibility of BEC systems with respect to changes of 
parameters permits, indeed, to investigate these phenom- 
ena in much an easier way than in their condensed mat- 
ter physics counterparts. On the other hand, BEC sys- 
tems are intrinsically nonlinear and the correspondence 
with usual condensed matter phenomena can be estab- 
lished mainly in the linear regime when the interatomic 
interactions are detuned to zero by means of external 
magnetic fields using the Feshbach resonance technique 
14j . The presence of nonlinearity, however, represents 



an additional resource for BEC systems which leads to 
interesting phenomena such as the existence of localized 
bound states which remain stable for infinite time due to 
a perfect balance between nonlinearity and dispersion. In 
presence of a periodic potential or a linear optical lattice 
(LOL), these states are known as gap-solitons (they have 
chemical potentials located inside the band-gaps of the 
underlying linear band structure), which can exist for 
both attractive and repulsive interactions 15Hl7j|. this 
last fact being only possible due to the presence of the 
periodic potential. In this context, it has been shown 
that gap-solitons are formed through the mechanism of 
modulational instability of the Bloch states, at the edges 
of the Brillouin zone of the underlying linear periodic 



system [l8(. Their existence in BECs with repulsive in- 
teractions was experimentally demonstrated in Ref. [l9l ]. 
Recently, the possibility of Bloch oscillations [j| , Rabi os- 
cillations [2(| and dynamical localization 0] of BEC gap 
solitons, in presence of constant or time dependent lin- 
ear forces (induced by accelerations of the OL), have also 
been considered. In these cases, besides the LOL, peri- 
odic spatial variations of the nonlinearity, also known as 
nonlinear optical lattice (NOL), have been used in order 
to avoid dynamical instabilities and to make the above 
phenomena long lived in the nonlinear regime. 

In higher dimensions, LOLs were shown to be effective 
in stabilizing localized states against collapse or decay, 
leading to the formation of stable multidimensional gap- 
solitons 
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[22j. In particular, it has been shown that 
a one-dimensional (ID) LOL in three-dimensional (3D) 
space does not allow to stabilize 3D gap-solitons. This 
soliton stabilization becomes only possible for attractive 
3D BEC subjected to the action of a 2D LOL [H, 
The possibility to stabilize solitons by means of NOL has 
been investigated mainly in the ID case, for which the 
mathematical properties have been studied in detail, for 
the ground state when considering solitons supported by 
ID NOL 241427 Las well as by combinations of linear and 
nonlinear OLs 28j. In the 2D case, it was shown in [29[ 
that, for conservative systems, a NOL in one direction 
by itself cannot give stable localized solutions in the case 
of attractive interactions. This fact remains true also for 
2D NOLs both for attractive and repulsive interactions. 

Since NOLs can be created by means of external mag- 
netic fields via the usual Feshbach resonance technique, 
or by optically induced Feshbach resonances (OFR) [30|, 
it is of interest to investigate whether crossed combined 
linear and nonlinear optical lattices (e.g., a LOL in one 
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direction plus a NOL in the orthogonal direction) allow to 
stabilize multidimensional solitons. This analysis can be 
relevant in the perspective of experimental observations 
of multidimensional solitons. In this respect, we remark 
that an experimental realization of a NOL was recently 
achieved (on the submicron scale) for a 174 Yb BEC us- 



ing the optical Feshbach resonance technique [31| . More- 
over, we recall that to date no multidimensional BEC 
solitons in LOLs have been observed. Besides, giving 
the possibility of exploring alternate methods for multi- 
dimensional soliton creation, crossed combined linear and 
nonlinear OLs can allow to extend to a 2D setting inter- 
esting transmission and reflection phenomena of soliton 
wavepackets, such as matter- wave optical limiting pro- 
cesses and bistability phenomena considered in the ID 
case for potential applications as matter-wave limiters, 
BEC mirrors or cavities and atomic switches 32 - 34 1 . 



The aim of the present paper is to investigate the ex- 
istence of 2D matter-wave solitons in crossed OLs con- 
sisting of a LOL in the x— direction and a NOL in the 
y— direction. In particular, we show that crossed linear 
and nonlinear OLs allow to stabilize 2D solitons both for 
attractive and repulsive interatomic interactions. This 
will be demonstrated by means of analytical and nu- 
merical approaches: analytically, with the variational 
method and a multi-Gaussian ansatz; numerically, with 
relaxation methods and direct numerical integrations of 
the Gross-Pitaevskii (CP) differential equation. Due 
to the lattice anisotropy, the solitons display elliptical 
cross sections, which in our variational approach (VA) 
are accounted by a multi-Gaussian ansatz with differ- 
ent parameters for each spatial directions. Existent VA 
curves and corresponding stability properties of multidi- 
mensional solitons are given for several parameter val- 
ues in terms of chemical potentials and total energies as 
functions of the number of particles, using the Vakhitov- 
Kolokolov (VK) criterion for the stability. The results 
are then compared with direct numerical integrations of 
the full GP equation. 

The paper is structured as follows. In section II, we in- 
troduce the model equations and discuss the physical im- 
plementation of a 2D crossed linear and nonlinear lattice 
by using spatial modulations of the scattering length. In 
Section III, we use a multi-Gaussian variational approach 
to derive our results, presented for chemical potentials 
and total energies as functions of the number of atoms, 
for the case of 2D solitons in 2D crossed OLs, for both 
attractive and repulsive interactions, where the stability 
is investigated by the well known VK criterion. In Sec- 
tion IV, the results of the VA are compared with the ones 
obtained by direct integrations of the GPE, using both 
relaxation in imaginary time and real time propagations 
to check the stability. Finally, in Section V, the main 
results of the paper arc briefly summarized. 



MODEL EQUATIONS 

Multidimensional BECs in 2D crossed linear and non- 
linear OLs are described in the mean field approximation 
by the following Gross-Pitaevskii equation (GPE): 



ih 



dt 



2m 



V|^-Acos(2fcx)^- 5 (y)M 2 V', (1) 



where i\) = ?/>(r, t) is normalized to the number of atoms, 
Vf, denotes the Laplacian in dimension T> = 2, m 
is the atomic mass, and Acos(2fcx) is a LOL in the 
x— direction, with strength A and lattice constant n/k. 
The function g{y) represents a NOL in the y— direction, 
produced either by spatially varying magnetic fields near 
a Feshbach resonance or by optically induced Feshbach 
resonances It the following we assume g(y) of the 

form 



g(y) = ,9o + .9i cos(2m/), 



(2) 



where go denotes the mean nonlinearity related to the 
mean s— wave scattering length a S Q, and g\ is the strength 
of a periodic modulation of the nonlinearity in the 
y— direction, having the period 7t/k. Note that the quasi- 
2D system is confined in the z— direction, with the ef- 
fective 2D scattering length being given in terms of the 
3D scattering length and some typical scale in the con- 
fined direction (3a, H(| . In this case, the parameters go 
and <?i are given in units of energy multiplied by some 
squared distance. The spatial modulation can be pro- 
duced by manipulating the scattering length with a laser 
field tuned near a photo association transition, e.g., close 
to the resonance of one of the bound p levels of the ex- 
cited molecules. Virtual radiative transitions of a pair of 
interacting atoms to this level can change the value and 
even reverse the sign of the scattering length. One can 
show that a periodic variation of the laser field intensity 
applied in the y— direction of the form I(y) = Iq cos 2 (k?/) 
produces a periodic variation of the atomic scattering 
length, such that a s {y) = a s o[l + aI/(5 + /)], where a s o 
is the scattering length in the absence of light, S is the 
frequency detuning of th e lig ht from the resonance, and a 
is a constant factor 2% |3jj. For weak intensities, when 
Iq 1 5 1, we have that the real part of the scattering 
length can be approximated by a s — a S Q + a s i cos 2 (fty), 
leading to a modulated nonlinearity of essentially the 
same form assumed in Eq. ([2]). 

It should be remarked, however, that the creation of 
a NOL in a BEC, by manipulation of the scattering 
length in this manner, also implies some spontaneous 
emission loss, which is inherent in the optical Feshbach 
resonance technique 37|-(39|. Such dissipative effects can 
be strongly reduced by using laser fields with sufficiently 
high intensity and detuned from the resonance It 
is also worth pointing out that, by using a laser field to 
control a magnetic FR the losses can also be essentially 
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reduced in comparison with optically induced FR. In 
particular, the experiment reported in (40| demonstrates 
that a laser light, near the resonance with a molecular 
bound-to-bound transition in 87 Rb, can be used to shift 
the value of the magnetic field where the FR occurs. By 
this way, it is possible to vary the scattering length on 
the optical wavelength scale, without having considerable 
losses (about two orders of magnitude lower then in op- 
tically induced FR experiments). We also remark that 
periodic variations of the scattering length, induced by 
the usual Feshbach resonance technique using spatially 
periodic external magnetic fields, would lead to similar 
conclusions. 

In the following, we adopt dimcnsionless units by 
rcscaling the space and time variables, such that the vari- 
ables in Eq. ([T]) are redefined according to t — > (h/E r )t, 
r — > r/k, with E r = h 2 k 2 /2m being the recoil en- 
ergy and r = {x,y, z}. The wavefunction is rescaled as 
ip(v,t) — > ku(r,t), in terms of which the Eq. ([T]) acquires 
the form 



. du 



dx 2 dy 2 



-Vi(x)u-r(y)\u\ 2 u, (3) 



where 



Vi{x) = ecos(2x), r(y) = x + 7Cos(Ay), (4) 



denote the linear and nonlinear OL interactions, respec- 
tively. In the above, e = A/E r , \ = gak 2 / E r = 8ira s ok, 
A = 2/i/fc, and 7 = gik 2 /E r , with gi assumed as a free 
parameter. 



VARIATIONAL ANALYSIS 

We consider localized solutions of Eq. ([3]) in the 
crossed linear and nonlinear OLs given by Eq. in 
the cases that the mean nonlinearity \ can be both neg- 
ative (repulsive interaction) or positive (attractive inter- 
action). In order to obtain analytical estimates for the 
soliton existence and stability, we use the variational ap- 
proach (VA) and look for solutions of the form 



u{x,y,t) = U(x,y)exp(-ifit), 



(5) 



where fi is the chemical potential and U{x, y) is a real 
function for the soliton profile. Due to the anisotropy 
of the crossed OL, having different natures and differ- 
ent strengths in the two directions, we expect elliptical 
cross sections for the soliton profiles. This feature can 
be accounted in our analysis of Eq. ([5]) by adopting the 
following ansatz for U(x,y): 



U(x,y) = Ae^ ax2+b y 2 ^ 2 , 



(0) 



where a and b are the parameters controlling the Gaus- 
sian widths in the two directions, with a/b the aspect 



ratio of the solution. Using this ansatz, the stationary 
GPE can be written as 

fiU + ^ + ^ + ecos(2x)U + [ X + ~f cos(\ y)]U 3 = 0, (7) 

from which the chemical potential and corresponding 
total energy are, respectively, given by 



y,N 



{\VU\ 2 - ecos(2x)U 2 - [ X + 7cos(A y)]U A }dxdy, 
2 x + 7cos(Ay) 



E = J j|W| 2 -ecos(2x)U 2 - 



-?7 4 j dx dy. 



The integrations in the above and following expressions 
cover the 2D phase space, from —00 to +00. From the 
ansatz ([5]) it follows that the normalized number of atoms 
is expressed in terms of the variational parameters A, a, 0, 
as: 



N 



U dx dy 



A 2 ir 



(8) 



The stationary GPE in ([7]) can be derived from the fol- 
lowing field Lagrangian 



L = - J {|VC/| 2 - [n + e cos(2x)]U 2 
— — [x + 7Cos(Ay)] U 4 \ dx dy 



(9) 



By substituting the ansatz (|6]) into the above Lagrangian 
and performing the integrations, in terms of the varia- 
tional parameters a, b and A, we obtain the following 
effective Lagrangian: 



L 



A 2 TT 



eff 



AVab 



A 2 



-2 



fi + ee 



2 

-l/a 



-A 2 /(86) 



with the corresponding total energy given by 



E — N 



a + b A 2 



-A 2 /(8&) 



ee 



-l/a 



(10) 



(11) 



Then, following from L e ff, by using the Euler-Lagrange 
equations for the parameters a, b, A, we can derive the re- 
lations for the chemical potential, total energy and num- 
ber of atoms. One can easily show that these relations 
can be written as 



b — 3a 



1 



N 

E 
N 



2 

4-7T a 2 - 2ee 



av ab 



X + le.' 



b-a ( 2\ 



(12) 
(13) 
(14) 
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with a, b solutions of the transcendental equation 



e at 



where 



b 7 46 2 



F(a) 



T = F(a), 



a 2 — 2ee 



(15) 



(16) 



Exact analytical solutions of Eq. (|15j) can be obtained in 
the following limiting. 

i) Case 7 = 0. In this case only the LOL in the x- 
direction is present and from Eq. (|15[) we obtain b = 
1/F(a). By substituting this value of b in Eqs. (|12|) and 
(|13[) . we get the following parametric equations, in the 
fi — N space, for the existent curves of a soliton: 



/' 

iV 



-a — e I 1 I e~ 

a , 

'l - e— e 

X 



(17) 
(18) 



Notice from Eq. fj 10[) . that this case becomes equivalent 
to the case with A = and 7 7^ 0, by replacing x m GU 
with x+7- The above equations coincide with the ones of 
the ID potential case investigated in Ref. [23| [compare 
the above equations with Eqs. (8) and (9) of this pa- 
per], where it was shown that the system can support 
2D localized solutions only for attractive interactions. 
The absence of a confining potential in the y— direction 
makes the solution to be uniform in this direction. But, 
interestingly enough, the solution is stable if the peri- 
odic boundary conditions are adopted in the y— direction 
(line-soliton). 

ii) Case x = 0- I n this case both the linear and the 
nonlinear OL are present but the mean nonlinearity is 
detuned to zero. Then, the Eq. (fT5")) reduces to 46 + A 2 = 
Ab 2 F(a) and b can be expressed in terms of the parameter 
a as 



1 + y/1 + \ 2 F(a) 
2F{a) 



(19) 



iii) Case e = 0. In this case only the NOL in the 
y-direction is present and from Eq. (|15p we have that 
F(a) = 1/a. The solutions for /x and N can be conve- 
niently expressed in terms of the parameter b as: 



a 
b 




X + 1? 



/' 

N 



2 V b 



-in 



X + je 



(20) 

(21) 
(22) 
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FIG. 1: (Color online) Variational approach results for 2D 
solitons, for attractive (upper panels, with \ = 0-5) an d re- 
pulsive (lower panels, with \ — —0.5) mean nonlinearity. The 
chemical potentials /j, are shown in the left panels, with the 
total energies E shown in the right panels, as functions of iV 
for increasing strengths of the linear OL e, as indicated inside 
the panels. Other parameter are fixed as A = 2, 7 = 1.0. All 
the quantities are in dimensionless units. 

In the absence of any confinement in the x— and 
y— directions (e.g., for e = and 7 = 0) the above equa- 
tions reproduce the known results of the 2D nonlinear 
Schrodinger solitons. In particular, for attractive inter- 
actions (x > 0), the soliton widths in the two directions 
become equal (a = b) and the number of atoms times 
the nonlinearity reduces to a constant, which coincides 
with the norm of 2D Townes solitons (4lj determining 
the critical threshold for collapse of a soliton in 2D as 
N cr = 4vr/x. 

In the general case the solution of Eq. (fT5|) must be 
found numerically. In Fig. Q] we show variational curves 
N, [i (left panels) obtained from the above equations for 
attractive (top) and repulsive (bottom) cases. Notice the 
non monotonic behavior with the appearance of a mini- 
mum (turning point) in N which gives a threshold in the 
number of atoms for the existence of the soliton. This 
is a characteristic of higher dimensional solitons of cubic 
nonlinear Schrodinger equations. It is associated to the 
so called delocalizing transition 42j; e.g., for values of N 
smaller than the critical value (minimum of the curve), 
the solution quickly decays in the uniform background 
and becomes fully delocalized. In the right panels of 
Fig. Q] we also show the corresponding plots in the E — N 
plane. Notice that, from the general definition, fj, = 
the chemical potentials correspond to the slopes of these 
curves at a given N. Also notice that the E — N curves 
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display cusps in correspondence of the turning points of 
the n~ N curves. 

The stability of the soliton solutions can be inferred 
from the observed plotted results b y m eans of the 
Vakhitov-Kolokolov necessary criterion |43[ , according to 
which stable solution always correspond to the branches 
for which dfi/dN < 0. We see that no solution can be 
stable for both attractive and repulsive interactions in 
presence of only the nonlinear OL (e.g. for e = 0). By 
increasing the strength of the linear OL, however, the 
slope of the curves changes from positive to negative. 
Our analytical study therefore predicts that stable soli- 
tons in crossed linear and nonlinear OL can exist not only 
for attractive interactions, a fact which is true also in ab- 
sence of the NOL in the y— direction as demonstrated in 
231 ] , but also for repulsive interactions (in absence of the 
NOL this last case would be impossible). It is also worth 
to note that the stable branches of the curve in the at- 
tractive case always lie well below the critical threshold 
for collapse, e.g. N < N cr , 

In the next section we show that these predictions are 
in good agreement with the results obtained by direct 
numerical integrations of the GPE. 



NUMERICAL RESULTS AND VA COMPARISON 

To check the above analytical predictions we have per- 
formed direct numerical integrations of the 2D GPE in 
Eq. (J3j> with crossed linear and nonlinear OL. The 
method we have used to find localized solutions is a 
standard relaxation algorithm in imaginary time prop- 
agation with a back-renormalization, which allows to fix 
the chemical potential /i and determine the correspond- 
ing number of atoms N 44 1 . 

Using this method we were able to follow these solu- 
tions in the parameter space even in the regions where 
they become unstable (hyperbolic states), which permit 
us to trace curves in the N — \i plane and to determine 
the existence of the thresholds in the number of atoms 
(turning points). Typical numerical curves for 2D soli- 
tons are depicted in Fig. [2j for both attractive (top pan- 
els) and repulsive (bottom panels) interactions for dif- 
ferent parameter values. Examples of soliton solutions 
determined with the relaxation method are shown in the 
left panels of Fig. |3l To facilitate the comparison with 
the VA predictions, the same set of parameters as the 
ones used in Fig. Q] are considered. We see that in the 
attractive case the agreement is very good not only quali- 
tatively but also from a quantitative point of view. Apart 
from the small shift of the curves, the VA correctly pre- 
dicts the existence of the turning point and the change 
of stability of the solutions. Notice that the agreement 
improves as the strength of the LOL is increased this be- 
ing a consequence of the fact that in a deeper LOL the 
soliton becomes more localized and better described by 
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FIG. 2: (Color online) Plots of 2D solitons in the (i - N 
(left panels) and E — N (right panels) planes, as obtained 
from numerical integrations of Eq. ©, for attractive (upper 
panels, % — 0-5) and repulsive (lower panels, \ = —0.5) mean 
nonlinearity and for increasing strengths e of the linear OL 
(see legends inside panels). Other parameters are fixed as in 
Fig. [T] Symbols on curves correspond to localized solutions 
selected for comparison with the VA density profiles in Fig. 
[3] and for stability checks by time propagation in Fig. [4] All 
the quantities are dimensionless. 



the Gaussian ansatz used in our variational analysis. For 
the repulsive case the agreement becomes more qualita- 
tive with a larger shift between the numerical and the 
VA curves, which is unaffected by the strength of the 
LOL. This discrepancy can be ascribed to the fact that 
for repulsive interactions the Gaussian ansatz used in the 
VA becomes less accurate because of the tunneling of 
the matter into adjacent potential wells (the condensate 
wavefunction cannot be localized into a single potential 
well and develops satellites in adjacent wells of the OL). 

In Fig. |3] we compare the sections of the density pro- 
files of the solitons corresponding to the points shown on 
the e = 1.0 curves in the corresponding top and bottom 
left panels of Fig. [2j The right panels refer to VA so- 
lutions for attractive (top) and repulsive (bottom) cases, 
while the corresponding profiles obtained from numerical 
imaginary time integrations of the GPE are shown in the 
left panels. We see that, for the chosen parameters, the 
agreement with our analytical predictions is quite rea- 
sonable. Also notice that for the attractive case (positive 
mean nonlinearity) the localization induced by the LOL 
in the £— direction is stronger than the one induced by 
the NOL in the y— direction. By changing the sign of the 
mean nonlinearity, however, keeping all the other param- 
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FIG. 3: (Color online) Sections of the soliton density taken 
at x = (dashed lines) and y = (continuous lines) for 
attractive x = 0.5 (top panels) and repulsive \ = —0.5 (bot- 
tom panels) interactions. Left panels refer to the numerical 
GPE imaginary time relaxation method calculations and cor- 
respond to the points shown on the e = 1.0 curves in the upper 
and lower left panels of Fig. at N = 6.02 {fx = -0.523) 
and N = 21.12 (y, = —0.96), respectively. Right panels show 
the soliton profiles predicted by the VA for the same chemical 
potentials as in the corresponding left panels. In this case, 
for fj, = -0.523, we have N = 6.2, a = 1.3 and b = 0.928 
(attractive case); and, for y = —0.96, we have N = 22.69, 
a = 1.93 and b = 2.6 (repulsive case). Other parameters are 
fixed as in Fig. [T] All the quantities are dimensionless. 



eters the same, the opposite situation occurs. 

The stability of these solutions has been checked by 
means of a real time propagation of the GPE using a 
Crank-Nicolson split-method and taking as initial con- 
dition slightly perturbed solitons. In the left panels 
of Fig. [2j we have selected four specific soliton solu- 
tions on branches of the curves, for e =2, with differ- 
ent slope close to the turning points, to check stabil- 
ity, e.g. (j, = — 0.7(diamond) and \i = — 0.55(star) in 
the attractive case(top frame); [i = —1.0 (sphere) and 
fi = — 0.55(trianglc) in the repulsive case (bottom frame). 
In Fig. [5] we show the dynamics obtained from numeri- 
cal integrations of the GPE with initial conditions cor- 
responding to these solitons with a small perturbation 
added. We see that while the solitons on the branch of 
the curves with negative slope are stable under time evo- 
lution, the ones on the branches with positive slope are 
unstable and quickly decay into an uniform background. 
These results are in agreement with what is expected 
from our VA analysis using the Vakhitov-Kolokolov cri- 
terion. From these studies we conclude that it is possi- 



FIG. 4: Time evolution of the x and y sections of the soliton 
density for attractive (top panels) and repulsive (bottom pan- 
els) interactions corresponding to the points on the stable (left 
panels) and unstable (right panels) branches of the e = 2.0 
curves in the top and bottom panels of Fig. [2] The chemical 
potentials of the stable and unstable solitons are y = —0.7 
and y = —0.55 for the attractive case, with y = —1.0 (left) 
and y = —0.55 for the repulsive case, respectively. In the left 
side of the frames we show the density scales, with the time 
evolution shown in the right side. Wider profiles correspond 
to y-sections. All the quantities are dimensionless. 

ble to create stable 2D BEC solitons using a ID linear 
OL in the x— direction and a periodic modulation of the 
scattering length in the orthogonal directions, both for 
attractive and repulsive interactions. 

An estimate of the parameters for possible experimen- 
tal observations of the above 2D solitons can be made 
by considering, for example, the case of a repulsive BEC 
with N ~ 10 5 atoms of 87 Rb loaded in a crossed com- 
bined OL consisting of a LOL in the X— direction with 
a period d = 2.2^1X1, generated by a laser field of wave- 
length n/k rj 0.830yum and strength A = 1.5 hx kHz 
(h =Planck constant), and a NOL in the y— direction 
created by an optically induced Feshbach resonance lead- 
ing to a spatial variation of the scattering length of the 
form a s (y) = a S Q + a si cos(2ny). The background scat- 
tering length, a S Q. is obtained by using an uniform 2D 
external magnetic field, with the value near the FR value 
Bfr = f007 G, which will give a s o ss 100ao, with ao 
denoting the Bohr radius. The modulation part of the 
scattering length a s i is induced by illuminating the sys- 
tem in the y— direction with a laser field of wavelength 
tt/k « 0.784/im, corresponding to the photo association 
wavelength of 87 Rb, to shift the value of a s i to values 
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Re(a s i/a s a) ~ ±(1 — 5)[40|. Such a crossed combined 
lattice is estimated to trap fundamental solitons (e.g soli- 
tons localized in a single potential well) containing up to 
few thousands of atoms. 



CONCLUSIONS 

By summarizing our present findings, we have shown 
that it is possible to stabilize two-dimensional BEC soli- 
tons by combining a linear OL, in one direction, with 
a nonlinear OL in the other orthogonal direction, where 
the NOL can be obtained from a periodic modulation 
of the scattering length. In particular, we have shown 
that in 2D crossed linear and nonlinear OLs families of 
2D solitons can exist and can be stable for both attrac- 
tive and repulsive interactions. We have determined ex- 
istent curves for chemical potentials and total energies 
in terms of the number of atoms, both by a Gaussian 
variational approach and by direct numerical integra- 
tions of the corresponding GPE. As shown, both varia- 
tional and full-numerical results have been found to be in 
good agreement. The stability of the solutions has been 
checked in both cases, by using the Vakhitov-Kolokolov 
criterion and by direct numerical real time integrations 
of the GPE. In this last case, for a few set of param- 
eters, the time evolutions of density profiles have been 
presented. These results open the possibility to observe 
two-dimensional localized matter waves in the presence of 
crossed linear and nonlinear optical lattices in real exper- 
iment, by using, for example, the technique considered in 
Ref. 0. 
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